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Abstract 

Starting from a general, microscopic fermion-to-boson mapping that preserves 

Hermitian conjugation, we discuss truncations of the boson Fock space basis. 

We give conditions under which the exact boson images of finite fermion 

operators are also finite (e.g., a 1+2-body fermion Hamiltonian is mapped to 

a 1+2-body boson Hamiltonian) in the truncated basis. For the most general 

case, where the image is not necessarily exactly finite, we discuss how to make 

practical and controlled approximations. 
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Pairwise correlations are often important in describing the physics of many-fermion sys- 
tems. The classic paradigm is the BCS theory of superconductivity where the wavefunc- 
tion is dominated by Cooper pairs which have electrons coupled up to zero linear momentum 
and spin; these boson-like pairs condense into a coherent wavefunction. Another example is 
the phenomenological Interacting Boson Model (IBM) for nuclei , where many states and 
transition amplitudes are successfully described using only s- and d- (angular momentum 
J = 0, 2) bosons, which represent coherent nucleon pairs. In both cases the large number of 
fermion degrees of freedom are well modelled by only a few boson degrees of freedom. On 
the other hand, however, despite some forays by Otsuka et al. 0, a rigorous microscopic 
basis for such phenomenological models is lacking. 

The basic problem is to represent the underlying fermion dynamics and statistics with 
a boson image amenable to approximation and numerical calculation. Considerable effort 
has gone into mapping fermion pairs into bosons [§]||. However these mappings typically 
suffer from a variety of defects. Most, such as the Belyaev-Zelevinskii || and Marumori |7j 
mappings give rise to boson Hamiltonians with infinite expansions, that is, /c-body terms 
with k — > oo. Convergence is slow even when "collective" fermion pairs are used. Finite 
but non-Hermitian boson Hamiltonians have also been derived, such as the well-known 
Dyson mapping [§J§. Such non-Hermitian mappings can mix 'physical' and non-physical or 



'spurious' states [T(| as discussed below. 

Less well-known are finite Hermitian alternatives to the Dyson mapping, obtained by 
mapping commutation relations of pairs of fermion creation and annihilation operators 



11| , [12| . These alternatives also have drawbacks: they either require a complicated, infinite 
operator (our "norm" operator defined below) to project out the correct Fermi statistics, 
or else the entire boson Fock space must be retained, which is far more unwieldy than the 
original fermion problem and which defeats the purpose of boson mappings. 

In this paper, by an alternate approach that maps matrix elements rather than commu- 
tation relations, we find that matrix elements of finite (specifically 1- or 2-body) fermion 
operators can be represented, in the full Fock space, by a finite boson image times an infinite 



boson "norm" or projection operator; because the image does not mix physical and spurious 
states projection is not needed. Furthermore the image of a Hermitian fermion operator is 



also Hermitian. This regains the result contained in p2| . We then give conditions under 
which one can obtain a finite and exact image in a truncated boson Fock space, and illus- 
trate with an example. In the most general case the exact image is not necessarily finite, but 
we show how to construct higher-order A;-body terms in a systematic and practical fashion 
which could lead to effective approximations. 

Consider a fermion space with 2Q single-particle states, and a fermion Hamiltonian H. 
The general problem is to solve the fermion eigenvalue equation 

H\%) = E P \%), (1) 

find transition amplitudes between eigenstates, and so on. To do this we require a many- 
body basis. Often the basis set for many-fermion wavefunctions are Slater determinants, 
antisymmeterized products of single-fermion wavefunctions which we can write as products 
of the Fock creation operators a], j = 1, • • •, 2f2 on the vacuum a| x • • • a] n |0) for n fermions. 
These states span the antisymmetric irreducible representation of the unitary group in 2Q 
dimensions, SU(20). But for an even number of fermions one can just as well construct 
states from N = n/2 fermion pair creation operators, 

N 



m = n A k i°) > ( 2 ) 

m=l 

with 

^^EK^H- (3) 

We always choose the Q(2Q — 1) matrices to be antisymmetric to preserve the underlying 
fermion statistics, and we choose for the normalization the trace tr (AqA^ = 5 a p. For this 
paper we represent generic one- and two-body operators by 

T = £%4a i; V = J2(^\VW)Al A ^ ( 4 ) 

ij [iv 



We begin with the straightforward mapping to boson states 



N 



hk>-too=n 0)> ( 5 ) 

m=l 

where the are boson creation operators. In conjunction with this simple mapping of states 
we construct boson representations that follow the philosophy of the Marumori expansion 
and preserve matrix elements of the fermion operators, for example introducing boson 
representations T B , V B , and most importantly H B such that (<f) a \ T B \<fip) = (ip a \ T 
{4> a I Vb = (ip a \ V \ ipp) , and (<p a \ H B \<pp) = (ip a \ H In addition, because wavefunc- 
tions of the form (0) do not form an orthonormal set, we construct the norm operator J\f B 
with the property (<f) a \N B 10/?) = (ip a \^p)- The construction of these operators, which will 



be given in detail in [13|, follow directly from the matrix elements which are found by gener- 



alizing the vector coherent state method of Reference [jTJJ. This procedure differs from the 
the usual Marumori expansion in that the latter does not have an explicit norm operator. 
With this mapping the fermion eigenvalue equation (|I]) becomes a generalized (because of 
the norm) boson eigvenvalue equation 

H B \%) = EpSfs \%) . (6) 

Because we have defined our boson operators so as to preserve matrix elements, the original 
energy spectrum of (0) is found in (H). However, because the boson space is much larger than 
the original fermion space, (||) also has additional spurious states that do not correspond to 
physical fermion states. These by construction have zero eigenvalue and do not mix with 
the physical space. 

When one constructs the norm operator |U| one finds it can be conveniently expressed 
in terms of the kth order Casimir operators of the unitary group SU(2f2), = 2 fc tr (P) fc , 
P = J2ar &i& T A ff A| (and so is both a matrix and a boson operator; the trace is over the 
matrix indices and not the boson Fock space) 

A^ = :exp(-£t^(7^: (7) 



where the colons ':' refer to normal-ordering of the boson operators. This result is also 
found in Expanding ([7]) in normal order one obtains the form [TJ[] 

oo k k 



Mb = i + E E < ( W ; M) n & L n ^h- 



(8) 



k=2 °"T m=l n=l 

The fc-body boson terms embody the fact that fermion pair creation and annihilation oper- 
ators do not have the same commutation relations as do boson operators, and act to enforce 
the Pauli principle. 

Although the representations T B ,V B are also complicated many-body operators similar 
in form to (H), we can write them in compact form, a result we have not seen previously in 
the literature, 



f B = 2E:tr [A CT TAtQ] b%Af B : 

<T,T 

Vb = E \ V \ E : { tr [ A « A l G ] tr [A„At G 

fl,U <7,T 

+4tr [A^AtPGA.AtG] } b%Af B : 



(9) 



(10) 



where G = (1 + 2P) _1 [13|. In general the boson representations given in (|7]), (|9|) and (|10"|) 
do not have good convergence properties, so that simple termination of the series in fc-body 
terms as in, e.g., (§) is impossible and use of the generalized eigenvalue equation (||), as 
written, is problematic. 

The explicit forms of <^), ( |I0"D suggests, however, that these representations factor in a 
simple way: T B = NbT b = T B Jv B and V B = M B V B = V B JV B , where the factored operators 
T B , V B , which we term the boson images of T, V, have simple forms. For example, a one- 
body fermion operator has a one-body boson image 



f B = 2Etr (A CT TAt) &t K , 
To prove factorization, one puts TbMb into normal order: 

2 E tr (A ff TAt ) {tfJf B b T - bl [Af B , K] } 



(11) 



(12) 



and uses [ 13 1 the completeness relation 



E ( A «L ( A «)/i' = | (^.i'^j' - ^'^.i') • ( 13 ) 

a 

and the resulting identities 

2 ^ tr(QA]jtr(A Q R) = tr(QR) - tr(Q T R), (14) 

a 

2 tr(QA^RA a ) = tr(Q)tr(R) - tr(Q T R). (15) 

a 

We can also show that [T B ,J\f B ] = from these identities. 

The two-body interaction V defined in (f|) can be rewritten in terms of products of two 
one-body fermion operators plus a remainder one-body by rearranging the fermion Fock 
operators. We can then map these one-body operators in terms of boson operators (10). 
Because these one-body operators commute with the norm operator, the interaction will as 
well. Normal ordering we get 



b l b » + 2 E E tr ( A - A l A -' A r A,4) b\b\,b T b T , 



(16) 



and in general one can find a image Hamiltonian H B = T B + V B . This result, and its relation 



to other mappings such as the non-Hermitian Dyson mapping, is found in Marshalek p^| . 

Thus any boson representation of a Hamiltonian factorizes: 7i B = J\f B H B . Since the 
norm operator is a function of the SU(2£7) Casimir operators it commutes with the boson 
images of fermion operators [|13|], and one can simultaneously diagonalize both H B and Mb- 
Then Eqn. becomes 

H B \%) = E' p \%). (17) 

where E' = E p for the physical states, but E' for the spurious states is no longer necessarily 
zero. The boson Hamiltonian H B is by construction Hermitian and, if one starts with at 
most only two-body interactions between fermions, has at most two-body boson interactions. 
All physical eigenstates of the original fermion Hamiltonian will have counterparts in (|17|). 
It should be clear that transition amplitudes between physical eigenstates will be preserved. 
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Spurious states will also exist but, since the norm operator Mb commutes with the boson 
image Hamiltonian Hb, the physical eigenstates and the spurious states will not admix. 
Identification of the spurious states is a serious though tractable problem, as Mb annihilates 
such states; furthermore spurious states can be shifted up in the spectrum through the use 
of the Park operator M =: (%: +4N(N — 1), which has zero eigenvalue for physical states 
and a positive definite spectrum for spurious states. 

A more critical question however is that of truncating the boson Fock space, by which we 
mean using states constructed from a restricted set of fermion pairs/bosons denoted by {a}; 



the operators in this space we denote by 



Mi 



B 



H 



to construct: for example, the norm operator is 



ft 



B 



B 



\k-l 



and so on, and are straightforward 



==«P (£M— tr PIT 



\k=2 



(18) 



where [P] T = J2&,f b\b? A^A\. These truncated representations still exactly preserve the 



matrix elements in the restricted fermion space: (0 a | Mb — (i>a i'p/ and so on. 

This is true even when the truncated set of fermion pairs represented do not form a closed 
subalgebra, a fact apparently overlooked previously |I0| . It should be evident that the 
truncated representations still do not mix physical and spurious states. 

Although the representations remain exact under truncation, the factorization into the 



image does not persist in general: 



H 



B 



H 



B 



Mb ■ This was recognized by Marshalek 



|12|1 and arises because the completeness relation fll3| ) in general is only satisfied if the 



complete Fock space is used. (An alternate formulation found in |l2j does not require 
the complete Fock space, but mixes physical and spurious states and so always requires a 
projection operator.) 

We can however find a sufficient condition such that a factorization 



n 



B 



H 



B 



A/" 



B 



(19) 



does exist, with Hb at most two-body; then with a further, stricter condition can guarantee 



Hb is Hermitian and commutes with 



Mi 



B 



First, consider a partition of the single fermion 



states labeled by i = (i a ,i c ), where the dimension of each subspace is 2Vt ai 2Vt c so that 
Q = 2Q a Q c . We denote the amplitudes for the truncated space as Ag and assume they can 
be factored, (A^- = (K% ja ® (At) icjc , with K^K = KK t = -L K T = (-1) P K, where 
p = (symmetric) or p = 1 (antisymmetric). 

Furthermore we assume the completeness relation (|T3|), which was crucial for proving 
that Hb = H B J\f B , is valid for the truncated space, 



£ 



Oj i'UA Al 



■3c 



(20) 



A necessary condition is that the the set of operators, Aa, A^ } form a closed 

subalgebra, as in the example given below. 

The norm operator in the truncated space then becomes 



B 



(21) 



where P = J2af b^bfA^Af so that [P] T = (^-) P- In this case the boson image of a 
one-body operator is the truncation of the boson image in the full space, 



B 



T, 



B 



T 



Mi 



B 



T, 



B 



2$> (A^TAt) b% . 



(22) 



(23) 



The representation of a two-body interaction can be factored into a boson image times 
the truncated norm, 



V, 



Vn 



Aft 



B 



T ' 



(24) 



but in the case where only (|20|) holds, however, Vb, while finite (1+2-body), is not simply 
related to [Vb] t as is the case for one-body operators and in fact is not necessarily Hermitian. 
We will discuss this general case elsewhere [^3|. Suppose one has the additional condition 



t " J jaj' c ,iai' c 



£</^H£(A,W ic (Aj 

W iaja 



(25) 



where the factor N a = Q a (2Q a + (— l) p ) is the number of pairs in the excluded subspace; 
while condition(^) looks complicated there are interactions that satisfy it, for example, two- 
body interactions constructed from one-body operators V = T & pT &l pi where T^p = A' a , As . 
When ( p5|) is satisfied then Vb is Hermitian and although Vb ^ [Vb]t they are simply related: 



Vb = J2(ct \V\ f) b\b f + 2f na £ (ji \V\ u) £ tr (A* A+ A^A} A.Aj,) b\b%b f b f , (26) 



with /n o = 4f^/iV a renormalizing the two-boson part of [Vb] t by a factor which ranges from 
unity (full space) to 2 for a very small subspace. 

The SO (8) and Sp(6) models [15] belong to a class of models which have a subspace 
for which (|20|) is valid. In these models the shell model orbitals have a definite angular 
momentum j and are partitioned into a pseudo orbital angular momentum k and pseudospin 
i, j = k + i. The amplitudes are then given as products of Clebsch-Gordon coefficients, 
(A\\ = + \ k rrii, krrij\K a M a ) (ifi^i /J>j\I a ^a), where K and I are the total pseudo 

orbital angular momentum and pseudospin respectively of the pair of nucleons. For the 
SO(8) model i = | and one considers the subspace of pairs with K — (p — 0), (A\)ij = 
( - 1+( -~ 1 ' > l (i fij\I a fi a ); in the Sp(6) model k = 1 and one considers the subspace with 
I — (p = 1), (Ag^ij = 1 (kmi,kmj\K a M a ). For these subspaces the completeness 

relation (|20|) holds, and the multipole-multipole interactions of [|15|] satisfy (|25|) and yield 
finite, Hermitian images as defined in (|26|) . 

Not all interactions satisfy (pB"D; for example, the pairing interaction never does [13 



except in the full space. In such cases the factored boson image is usually non-Hermitian 
|T3|. Since the norm operator and the boson representation of Hermitian fermion operators 



are Hermitian, and under factorization 



Hb 


T = H B 


'Mb 


T 


Hb 


t 

T 


'Mb 



T H B , the 



boson image in the truncated space is Hermitian if and only if it commutes with the norm. 
If the image commutes with the norm then they can be simultaneously diagonalized and the 
image does not mix physical and spurious states. If the image does not commute, however, 
as is the case for the non-Hermitian Dyson mapping, then it can mix physical and spurious 
states. This is because the norm is a singular operator with null eigenvalues; even though the 
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representations have zero matrix elements between physical and spurious states, by factoring 
out the norm one essentially is dividing zero by zero which can lead to mixing. 

For truncation schemes where the conditions ( p0| ) or ( p5| ) do not hold, it is better to 
define a hermitian image Hamiltonian 



-1/2 



-1/2 



Mb „ [H B ] T Mb „ ■ (27) 



T 

1/2 . 



Because the norm is a singular operator it cannot be inverted. Instead Mb is calculated 
from the norm only in the physical subspace, with the zero eigenvalues which annihilate the 
spurious states retained. Then Hb does not mix physical and spurious states. When the 
image commutes with the norm this is equivalent to our previous definition. 

Hb is not necessarily finite; it may have an infinite expansion which for purposes of 
discussion we sketch as 

H T = 0j&t& + 6 2 b ] b ] bb + 9zb ] b ] b ] bbb + ... (28) 
We now argue that calculation of the coefficients 6k for the fc-body term can be done in a 

r ~ 1 -1/2 

straightforward manner. Because the norm operator is of the form (]^), the operator Mb t 

must be of a similar form. Further the fc-body term can be calculated by diagonalizing [Mb] t 

r ~ r i - 1 / 2 

in the truncated fc-body Fock space, and 6k at most requires the fc-body terms of Mb 
and \Hb\t- Then one can consider the convergence of the series ( p8|) as a function of k. A 
rough estimate is that, for an iV-boson Fock space, one can truncate to the £-body terms if 
for k > £, 8k is sufficiently small compared to 8e x (N — k)\/(N — £)\; the strictest condition 
is to require ^ < ^/(fc-l)!. 



Under a finite termination of ( 28[) Hb can again mix physical and spurious states, but if 
the remainder is small the mixing may be small. (If the truncation in the boson Fock space 
is severe enough the boson Fock space may not even contain any spurious states, as for the 
SO (8) model below half-filling.) 

We emphasize the utility of this "linked-cluster" expansion similar to that of ]TB|] . Sup- 
pose the series (^) converges fast enough that truncation to 4-body terms is sufficient even 

10 



for states with 5 or 6 bosons by the criterion given above. If one takes just J = 0, 2 bosons 
(as for the interacting boson model), by partitioning via angular momentum the largest 
matrix that needs to be inverted to calculate the 4-body terms is of dimension 4. This ap- 
proximate image Hamiltonian, although calculated in a 4-boson system, would be applicable 
for iV-boson systems with N ^> 4, failing only when a significant number of spurious states 
enter the Fock space. 

As a final example, consider the pairing interaction in the seniority-0 subspace, where 
there is only one state, (S^) N |0), with = (4fi)~ 1 / 2 J2i e^ajal; here % is the time-reverse 
of the state i, 6{ = ±1 and e% = — e*. Under naive truncation the pairing interaction, 
V p = —GS^S, gets mapped to 

~V P 

which does not give the correct ground state energy. Because condition fl2"5|) does not hold 
one cannot use ( pip to obtain the image; instead (^) yields the correct answer: 

V£ = -G (sU - ^shs 

Applications to realistic systems will be discussed in future investigations. 
This work was supported by the U. S. Department of Energy. 
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(29) 
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